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Abstract 



We calculate, in the framework of multi Higgs doublet models, the CP violat- 
ing transverse tau polarization, P±, in inclusive semileptonic B meson decays. 
We find that P± diverges at 0(l/m%) of the heavy quark expansion. We dis- 
cuss the physical reasons for this divergence and show how to regularize it. 
We find large corrections that can suppress the CP asymmetry by as 

much as 30%. The maximal allowed asymmetry is acp ~ 0.34. We discuss 
how the allowed range is expected to change in the future. 
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I. INTRODUCTION 



Multi Higgs doublet models (MHDM) allow for new CP violating phases in charged 
scalar exchange ffl. In models with natural flavor conservation (NFC) |J at least three 
Higgs doublets are needed for the scalar sector to allow CP violation The effect cannot 
be large in CP asymmetries in neutral B decays or in CP violation in the K — K system [|J. 
However, there could be significant contributions to processes for which the SM prediction 
is small: the electric dipole moment of the neutron H, CP asymmetries in top decays ||, 
and CP violation in semileptonic decays. 

As triple-vector correlation is odd under time-reversal, the experimental observation of 
such correlation would signal T and - assuming CPT symmetry - CP violation. In three 
body semileptonic decays the lepton transverse polarization is such an observable. It cannot 
be generated by vector or axial- vector interactions only 0, so it is particularly suited for 
searching for CP violating MHDM contributions. The muon transverse polarization has 
been studied in K — > tt fi v decays |7|,|8j . Since in models with NFC scalars couple more 
strongly to heavier fermions, the expected signals are stronger in heavy quark decays |9|-|T2"||. 

For top decays the quark level description is valid, as the top quark is expected to 
decay before it hadronizes. This is, however, not the case for B decays. Recently, it has 
been observed that inclusive semileptonic decays of hadrons containing a single heavy quark 



allow for a systematic, QCD-based expansion in powers of 1/rriQ [pL3 1 . The heavy quark 
limit coincides with the free quark decay model and there are no corrections to this result 
at order 1/mg Rl3,14|. The leading nonperturbative corrections are of order I/ttlq and 
depend on only two hadronic parameters, which parameterize forward matrix elements of 
local dimension-five operators in the heavy quark effective theory (HQET). These corrections 
have been computed for a number of processes [|14"|--|18|j . 

Using the combination of the operator product expansion and HQET, semileptonic B 
meson decays into a tau lepton have been studied in the framework of the SM fll6l , |17|l and 
MHDM [[L8|]. In this paper we use this method to study the tau transverse polarization 
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in the framework of MHDM. In addition to finding the correction to the spectator model 
result, this calculation has some interesting technical aspects. The nonanaliticity of the 
triple product correlation on the boundaries of the Dalitz plot gives rise to infinities in the 
calculation. Such divergencies do not appear in the calculations of either cross-sections or 
longitudinal polarization. We discuss their physical meaning and how to eliminate them. 



II. THE CALCULATION 



We consider a general multi Higgs doublet model with NFC (for a recent analysis and 



notation see [fOJ)- In order to have observable CP violation, the lightest charged scalar has 



to be much lighter than the heaviest one P0|. Here we assume that the all but the lightest 



charged scalars effectively decouple from the fermions. The Yukawa couplings of the lightest 
charged scalar to up-type quarks, down-type quarks and charged leptons are determined by 
the parameters X, Y and Z, respectively. The terms in the effective Lagrangian relevant 
for B — > r v X decays are 

AG F 



L = -V r 



cb ■ 



V2 
where 



l^P L b){f ltl P L p T )-R{cP R b){fP L v T ) , (2.1) 



R = r 1 m T m Y hl r 2 = XZ*/m 2 H , (2.2) 

Pr,l — o(l i 75), and m# denotes the mass of the lightest charged scalar particle. The 
running quark masses are denoted by an upper index Y. The first term in Eq. ( [2.1D gives 
the SM contribution, while the second one is that of the charged scalar. We neglect a term 
proportional to Y m%: first, it is suppressed by the mass ratio /m%; second, Y is bounded 
to be of 0(1), while X can be large |L9] . 



The tau transverse polarization, P±, in the decay B — > tvX is defined as the tau 
polarization component along the normal vector of the decay plane. It is given by 

P± = *r^M , (2.3) 
\Pt x Px\ 

.3 



where s T is the tau spin three vector and p T (px) is the three-momentum of the tau lepton 
(hadron). We define T + (T~) as the rate of finding the tau spin parallel (anti-parallel) to 
the normal vector of the decay plane. The total rate r = T + + T~ has been calculated in 
fl"8|| . We then have to calculate only r + — r~. 

The techniques of the calculation are described in detail elsewhere [|15||, so we give only 



the necessary details, followed by the results of the computation. The asymmetry arises 
only from the interference between the W mediated and the Higgs mediated diagrams. The 
interference term in the hadronic current can be decomposed as 

(2tt) 3 £ 5\p B -q-px) (B v \ J ] \X) (X\ J» \B V ) =v»U x + q» U 2 , (2.4) 
x 

where 

j» = cYP L b, J = —cmlXP R b J (2.5) 

m H 

g M is the momentum of the lepton pair, and is the four-velocity of the decaying B meson. 
We work in the B rest frame, where 9 is the angle between p T and px, and the tau spin 
projection is along the four-vector s = (0, 0, sgn[sin#], 0). Upon contraction of the lepton 
current with W^, the term proportional to U2 vanishes and the U\ part yields the differential 
asymmetry 



- — z—^- = ^ cb \ — ^ mb Ui ImR J xy(q 2 — ml) — x 2 mi — (q 2 — ml) 2 . 
aq z ax ay Ion* ' 



(2.6) 



where we defined dimensionless variables 



2E T 2E U ^ q 2 m c A m T 

V = , x = , q = — m c = — , m T = , (2.7) 

mb mb mi mb mb 



and E T (E u ) is the energy of the tau (neutrino). The OPE yields for U\ |TL8 



Ux = 5(z) + {[5(Ai + 3A 2 )(x +y) - 36A 2 ] S'(z) + X^q 2 - (x + yf] 5"{z)} , (2.8) 



where Ai and A 2 are parameters of the HQET, z = 1 + q 2 — m 2 — y — x, and the prime means 
d/dz. 



In order to obtain the total asymmetry we have to integrate ( |2.6|) over x, y, and q 2 . This 
integral diverges because the square root function is not analytic. The derivatives resulting 
from the 5' and 5" terms yield negative powers of the square root in ( |2.6| ) to be evaluated 
precisely at their pole. Physically, the reason for this non-analyticity is that the orientation 
of the decay plane changes when the tau and the hadron momenta become collinear. At 
these points in phase-space (corresponding to the boundary of the Dalitz plot in the y — q 2 
plane) the decay plane is not well-defined, so neither is the transverse polarization. In order 
to regularize this divergence, we have to rely on the physical picture. Experimentally the 
region near the boundary of the Dalitz plot is not accessible: any finite resolution gives a 
region in which the transverse tau polarization is not measurable as even the decay plane 
cannot be reconstructed. This suggests to regularize the divergencies by a simple cutoff, 
and then take the limit of the cutoff approaching the border of the Dalitz plot. Another 
possibility is to replace the square root by an analytic function, which means physically to 
make the flip of the orientation of the decay plane smooth. It is convenient to apply the 
following replacement 



>//(*, V, <f ) - y/f(x,V,q 2 ) [l - e~ nf } J , (2.9) 

and interchange the order of the integration with the n — > oo limit. This procedure elim- 
inates the infinities. In fact, this prescription is only sensitive to the derivatives of the 
regularized function at the points where the original square root function vanished (these 
points constitute the border of the Dalitz plot). Since at these points all derivatives of 
Q2.9| ) vanish, this should reproduce the cutoff result. We checked that in the m c — > limit 
these two procedures indeed give identical results. (The calculation using a cutoff becomes 
unmanageable when we keep m c finite.) Technically the replacement ( |2.9| ) is much simpler 
than doing the full calculation with a cutoff. 

After the integration we obtain for the CP violating asymmetry 



where 



•57T 



r_L = — { c A 

35 I 



6ml 



A 



2m 2 . 



2 ^3 



1 + 4m c — Arh c — rh c + 4m T (l + 5m c + m c ) — 4m r (l — m c ) — m T 
3(1 — m c ) 2 (8 + 5m c + m 2 ) + 9m T (8 — 3rh c — 4m 2 — ml) 
+4m 2 (8 - 9m c - 6m 2 ) + 9m^(l - m c ) + 3m^ 
3(— 16 - 13m c + 9m 2 + 15m 3 c + 5m 4 ) - 9m T (16 + 29m c + 20m 2 + 5m^) 
-4m 2 (2 + 45m c + 30m 2 ) + 45^(1 - m c ) + 15m*] } , (2.11) 



and C = (1 — m c — m T ). In the limit Ai, A2 — > 0, corresponding to the free quark decay 
model, our results agree with Ref. Jll| . (Although Ref. [[II] works in the tau rest frame, our 
definitions of the asymmetry are identical.) 



III. DISCUSSION 



In our numerical analysis we use the same input parameters as in Ref. [18]. Let us further 
define: 

r . -^tt J\Re(XZ*)\ r . J\\m(XZ*)\ , , 

r Re = ^/ Re r 2 =^-! ^ ^ , r Im = ^ Im (r 2 )| = ^ i <± . (3.1) 

m H m H 

The quantity ri m enters the asymmetry acp through its numerator quadratically, and the 
denominator in a more complicated manner. The dependence of acp on rR e is only through 
the total decay rate. In our analysis we use only the allowed values of r Im and r Re , and we 
discuss the present constraints on them below. By varying these parameters within their 
allowed ranges, we find the allowed values of acp as a function of r\ m . Our result is given 
by the shaded region between the solid lines in Fig. 1. 
A few points are in order regarding this result: 

a. The 1/m 2 corrections to the free quark decay model decrease acp for all values of 
ri m . The numerator is suppressed by about 25%. For small values or ri m the denominator is 
also suppressed and the total correction is small. However, for large values of r\ m the total 



branching ratio is enhanced [IS] and acp can be suppressed by as much as 35%. 



b. The uncertainty in our result (corresponding to the width of the shaded region in 



Fig. 1) takes into account the theoretical uncertainties as in Ref. ||18|| . 

c. The improvement of our calculation over the spectator model originates from large 
1 jm\ corrections, and from using rrib and m c as determined by the HQET rather than 
treating them as independent input parameters. To illustrate this, we plotted with dashed 
lines in Fig. 1 the prediction of the free quark decay model corresponding to 1.4 < m c < 
1.5 GeV and 4.6 < m b < 5 GeV. 

d. Even in the SM "fake" asymmetries arise as a background due to CP conserving 
unitary phases from final state interaction (FSI) and from real intermediate states. These 



are expected to be less than 1% [^TJ. Furthermore, they can be removed using the fact 



that asymmetries due to CP violating phases have different sign in B and B decays, while 



unitary phases yield asymmetries with the same sign |22]. Therefore, in our analysis we 
ignore them. 

The strongest bound on rj m comes at present from the measurement of the (unpolarized) 
B — > t v X branching ratio. The total decay rate T has the form 

r = W-IReR + H \R\ 2 (3.2) 

where W (I, H) denotes the SM (interference, Higgs) term. This is a circle in the complex 
{ r R e , r im} plane, centered on the r Re axis. In our analysis we set the value of r Re to allow 



for maximal 7i m subject to the experimental constraint. From the measurements [23-25 



BR (B ^tvX) = 2.98 ±0.47%, 

BR(B ^ivX) = 10.43 ±0.24%, (3.3) 



and using Ref. JT^] we obtain at the 2a level 

rim < 0.47 GeV- 1 . (3.4) 



Due to the negative sign of the interference term in ( |3.2| ), the strongest possible bound 
from the total branching ratio that can be obtained in the future is about 7i m < 0.35 GeV -1 . 



The longitudinal polarization of the tau in the same decay has a better sensitivity [|T8|j26 
down to about rj m < 0.25 GeV -1 . Note, that from just one CP conserving observable no 
determination of ri m is possible, since deviation from the SM prediction can also be induced 
by rR e . In principle, however, by comparing two observables with different dependence on 
r Re and r Im one may be able to detect a nonzero r Im . 



The purely leptonic decay B — > r v is also sensitive to r\ m . Using Ref. |27j] and the recent 
ALEPH result [E3I 



BR (5 -> rv) < 1.5 x 10" 



(90% CL) 



together with the f B > 140 MeV and \V ub \ > 0.0024 bounds, we get 

r Im < 0.58 GeV" 1 . 



(3.5) 



(3.6) 



Due to the dependence on ri m and the large theoretical uncertainties in the calculation we 
estimate that this bound cannot become stronger in the future than r\ m < 0.3 GeV -1 . 

The bound on ri m from P± of the muon in K — > it fi v decays is obtained by converting 



the measured quantity, Im£, to rj m [g8 



|Im£| 



2 2 
"im m K 



(3.7) 



The theoretical uncertainty in this relation enters through the calculation of certain form 
factors. Im £ has been measured in both charged and neutral kaon decays. The background 
of these measurements due to CP conserving FSI phases are different. For the charged kaon 
decay K + — > 7r° fi + u, FSI are second order in the electromagnetic interaction and therefore 
very small: (9(10~ 6 ) In neutral kaon decay Kl — > 7r~ /i + u, there are two charged 

particles in the final state, and so FSI is first order in the electromagnetic interaction and 
can contribute to Im£ at the 0(1O~ 3 ) level. Thus, in order to set a bound on r Im we use 
7T _ n + v) w lm£ FSI +rf m m 2 K instead of Eq. Q, where lm£ FSI w 0.008 p|,|0 



Im£(K z 



The PDG data 25 



Im£(K L -> 7T" /x + i/) = -0.007 ± 0.026 . 
Im£( J fs: + -> 7T° /x + v) = -0.017 ± 0.025 . 



(3.8) 
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yields at the 2a level 

r Im < 0.47 GeV" 1 . (3.9) 

Numerically, the bounds ( p.4[ ) and ( |3.9| ) are the same. However, we consider the bound ( p.9| ) 
less reliable as it is sensitive to the uncertainties in the form factors and in the FSI phases 
in Kl decay. 

In the near future the E246 experiment at KEK will measure the transverse muon po- 
larization in K + — > 7t° ji + v decay. They are aiming at a sensitivity to P± of about 5 x 1CT 4 
which corresponds to ri m ~ 0.08 GeV -1 [p8fl . From Fig. 1 we see that in order to achieve a 
better sensitivity to rj m from the measurement of P± of the tau in B — > r v X decay, one 
needs to measure P± at an accuracy of better than 2%. It is important whether P± of the 
tau can be reconstructed at all from its decay products in B decays, either at LEP, at B 
factories, or in hadron colliders. To our knowledge, this issue is not yet settled and further 
experimental studies are called for. 

One may naively expect that for the same CP violating parameters, the tau polarization 
in B decays would be mf ) m T /m s m^ ~ 500 times larger than the muon polarization in K 
decays. However, our calculation gives only about a factor of 40 enhancement. The scalar 
matrix element in exclusive K decays is enhanced by a factor of (mic/m s ) 2 ~ 7.6. The 
finite mass effects in the phase space integrals reduce a>cp by about 45% in B decays, while 
by only 30% in K decays. Finally, the l/m^ corrections reduce acp in B decays by up to 
35%. In addition, the experimental measurement of the transverse muon polarization in K 
decays is simpler than measuring the tau polarization in B decays. Thus we conclude that 
although most of the MHDM parameters are easier to probe in processes involving heavy 
fermions, it may be easier to constrain ri m in kaon decays. 

It is interesting to compare our results to Ref. |12|], which considers exclusive B decays. 
There it was found that the allowed ranges for acp extend from 0(35%) in B — > D*tv 
decay up to 0(1) in the B — > D r v channel. We believe that the suppression that we found 
in the inclusive channel is likely to apply to the exclusive ones as well. On the one hand, it 
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would be easier to measure acp in inclusive decays (larger statistics, no need to reconstruct 
specific hadronic final states) than in exclusive channels, and the exclusive calculation also 



suffers from larger theoretical uncertainties. On the other hand, the results of Ref. |12 
suggest that the asymmetry in the B — > D t v decay may be enhanced. 

Finally, our analysis holds with minor modifications for inclusive B — > /iuX and D — > 
fi v X decays as well. However, the possible signals are smaller and harder to measure 
experimentally than in the B — > t v X and K — > n /i + 1/ decay modes discussed above. 



IV. SUMMARY 

We studied the transverse tau polarization asymmetry in inclusive B —> t v X decays 
in MHDM. We use the heavy quark expansion to incorporate the nonperturbative I /ml 
corrections. Unlike previous cases, in the present calculation divergencies appear that have 
to be carefully removed. The l/m^ corrections turn out to be important and reduce the 
asymmetry by as much as 30% in certain regions of the parameter space. We studied the 
present and possible future experimental bounds on the CP violating parameter of the 
Higgs sector, \\m(XZ*)\/m 2 H . The strongest bound comes at present from the measured 
B — > tv X branching ratio. In the future - even if a new source of CP violation will not 
be discovered by E246 at KEK - the bound on the measurement of the muon transverse 
polarization in K — > tt \i v will provide the strongest constraint. 
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FIGURES 




FIG. 1. acp as a function of n m = y^Im (r 2 )| = -\/|Im (XZ*)\/rriH- The shaded area between 
the solid lines is our result. The area between the dashed lines gives the free quark decay model 
result. 
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